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The Kinematic and Chemical Properties of a Potential Core-Forming Clump: 

Perseus Bl-E 

S. I. Sadavoy^, Y. Shirley^’^, J. Di Francesco^, Th. Henning^, M. J. Currie'^, Ph. Andre^, S. 

Pezzuto® 

ABSTRACT 

We present and (1-0), (2-1), and (3-2) maps towards the core-forming 

Perseus Bl-E clump using observations from the James Clerk Maxwell Telescope 
(JCMT), Submillimeter Telescope (SMT) of the Arizona Radio Observatory, and IRAM 
30 m telescope. We find that the and line emission both have very complex 

velocity structures, indicative of multiple velocity components within the ambient gas. 

The (1-0) transitions reveal a radial velocity gradient across Bl-E of ~ 1 km s“^ pc“^ 
that increases from north-west to south-east, whereas the majority of the Perseus cloud 
has a radial velocity gradient increasing from south-west to north-east. In contrast, 
we see no evidence of a velocity gradient associated with the denser Rersc/iel-identified 
substructures in Bl-E. Additionally, the denser substructures have much lower systemic 
motions than the ambient clump material, which indicates that they are likely decou¬ 
pled from the large-scale gas. Nevertheless, these substructures themselves have broad 
line widths (~ 0.4 km s“^) similar to that of the gas in the clump, which suggests 
they inherited their kinematic properties from the larger-scale, moderately dense gas. 
Einally, we find evidence of depletion only toward one substructure, B1-E2, which 
is also the only object with narrow (transonic) line widths. We suggest that as prestellar 
cores form, their chemical and kinematic properties are linked in evolution, such that 
these objects must first dissipate their turbulence before they deplete in CO. 


1. Introduction 


Stars form within prestellar cores, which are small (< 0.1 pc), cold (< 10 K). and dense (> 
cm“^) objects inside molecular clouds (e.g., IWilliams et al.l 120011 : IPi Erancesco et al.M2nn7l ). 
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addition, prestellar cores are dynamically q uiescent both in their internal velocit y dispersion s 


Rosolowskv et al.l 120081 : iPineda et alJ 12010 1 and in their systemic motions (e.g., iKirk et al.ll201 


and most prestellar cores have similar centroid velocities between their densest, central mate rial and 
their moderately dense envelope (e.g., Walsh et al.ll2004 : Kirk et al. 2007 : Andre et aLlboOT l. These 
observations suggest that ( 1 ) prestellar cores have homogenized their relative motions compared to 
the surrounding, ambient cloud and ( 2 ) the densest gas material within these cores are coupled to 
(e.g., move quiescently with) their surrounding envelope. 

Prestellar cores also have significant chemical differentiation, where the higher density environ¬ 
ments are best probed with nitrogen-bearing molecules or deuterated tri-hydroge n molecules and the 


lower density envelopes are best traced with carbon-bearing molecules (e.g., see lDi Francesco et al 
20071 ). This differentiation can be attributed to different effective excitation densities between 


molecules (IShirlevll2015l ) and to different temperatures for molecular “freeze-out”, where certain 
molecular species adsorb onto dust grains and deplete from the gas ph ase. In particular, carbon- 


bearing molecules adsorb easily in co ld, dense prestellar cores (e.g., iBergin et al.l Il995l . 


2002 


Kramer et al.lllOOOl : ICaselli et al.l 119991 ) , and observation s of such objects often show absorption 
at 4.67 /im attributed to CO-ices (e.g., IChiar et al.lll994l ). Additionally, the formation of carbon- 
rich icy mantles on dust grains provide s a surface for larger, organic molecules to form via the 
irradiation of those icy mantles (e.g., see lHerbst &: van DishoeckI 120091 1 . 


Therefore, prestellar cores have distinct dynamical and chemical signatures compared with 
their parent cloud. Most clouds are actively forming stars or contain a current generation of dense 
cores, and thus, can no longer be used to constrain the initial conditions for their formation. As 


ditions, such as radiation feedback (e.g., ^ 

3at 

2009b: ( 

Dffner et al. 2009: Krumholz et al. 20101. 

magnetic fields (e.g.. Li & Na 

camura 

200 

1 

3asu et al. 

2009a 1 

)), different turbulent properties 

(e.g.. Klessen et al. 

2000: Bate 

2009a: 

Heitsch et al. 20111. or a combination of these effects (e.g.. 

Price Sz Bate 2009: 

Commercon et al. 

2011 : F 

’adoan &: Nordlunc 

1 20111. Thus, observations of a 


young core-forming region are key to constraining better the physical processes that cause cores to 
condense out of the larger-scale ambient cloud material. 

The Bl-E clump in the Perseus molecular cloud is an ideal laboratory to study core for¬ 
mation in relative isolation. Bl- E is the sole den se {Ay > 5 mag) region i n Perseus without 


young stellar objects (Y SOs; e.g.. lEvans et al.ll2009li and pre stellar cores (e . g., 


Enoch et al. 


2006 


Di Erancesco et al.l 120081 ). Nevertheless, using Herschel data. ISadavov et al.l (1201211 identified nine 
distinct substructures in Bl-E. Eollow-up NH 3 (1,1) observations of these substructures showed 
that most of them had supersonic non-thermal motions with only one source, B1-E2, having more 
quiescent, transonic non-thermal motions. Since the NH 3 emission from these substructures wa s 
weaker than the typical line peaks associated with prestellar cores (e.g., iRosolowskv et al.ll2008l ). 
the Bl-E substructures appear to be in an unique state where their dense gas is both minimal 
and turbulent. Thus, the Bl-E clump is an excellent candidate for a core-forming region that is 
fragmenting for the first time, such that its Herschel-identiHed substructures represent early-stage 
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core precursors. 

For this paper, we mapped the entire Bl-E ciump in and emission of their J = 1—0, 
J = 2 — 1, and J = 3 — 2 transitions. We use these data to characterize the kinematic and chemicai 
properties of the Bl-E ciump and its substructures. This paper is outiined as foiiows: in Section 
[21 we present our and observations; in Section |3l we describe the kinematic properties 

of Bl-E and its substructures using our hne observations; in Section we use radiative transfer 
modeiing to determine the chemicai properties across Bl-E; in Section [5l we discuss the observed 
kinematic and chemicai signatures, and we propose an evoiutionary picture for core formation; and 
hnaiiy, in Section [6l we present our conciusions and summarize the paper. 


2. Data 


2.1. Herschel 


The Perseus moiecuiar cloud w as observed by Herschel as part of the Herschel Gould Belt 
SurvejQ (HGBS; Andre et al. 2010 1. with the Western field mapped in Eebruary 2010 and the 
Eastern field mapped in February 2011. Both Perseus fields wer e observed using the PACS/SPIRE 
paral lel mode scanning technique with a 60 arcsec s“^ scan rate ( Sadavov et al.ll2012l : jPezzuto et al 


2OI2I : Pezzuto et al. in preparation). The data were reduced using HIPE version 7 and modified 
standard reduction s cripts by M. Sauvage (PAGS) and P. Panuzz o (SPIRE) in the same manner as 
Sadavov et al.l (|2014l l. We used scanamorphos 11 (|Rousselll2013l l to produce the final maps. 


We corrected the zero-point of fset at each waveleng th for both the Eastern and Western fields 
using Planck observations following [Bernard et al.l (I2OIOI ). Since both fields maps overlapped toward 
Bl-E, we mosaicked the two helds using the Starlink task wcsmosaic and bilinear interpoloation. 
Eor the 160 — 500 fiui data, the two helds had excellent agreement in their mosaics. Finally, we 
convolved the mosaics to a common resolution of 36.3" (the 500 /rm beam) and a common grid of 
14" pixels. 


2.2. James Clerk Maxwell Telescope 

We observed the Bl-E clump with the Heterodyne Array Receiver Programme (HARP) de¬ 
tector at the James Glerk Maxw ell Telescope (JCMT) in the same manner as the JGMT Gould 
Belt Survey ([Buckle et al.l l2009l ) . The HARP detector was tuned to map Bl-E simultaneously 


m 


I3r 


CO (3-2) and C^°0 (3-2) hne emission at 330.58796 GHz and 329.330545 GHz, respectively. 
Observations were taken using position-switching, orthogonal raster scans with 1/4 array spacing 
over a 25' x 14' region. The backend was conhgured to 4096 channels over 250 MHz, resulting in 


^http: //gouldbelt-herschel.cea.fr/ 
























- 4 - 


~ 0.055 km s“^ velocity resolution. Bl-E was observed on 13 July 2011 and 12 January 2012 for 
~ 2.2 hours of observing time (~ 20 % of the total project time). We selected the coordinates 
03:35:47.7, 31:45:30.5 (J2000.0) for our off-position. 


The HARP dat a were processed using the Auto-Correlatio n Spectral Imaging Sy stem (ACSIS; 


Jenness et al.l l2008l ) an d reduc e d usi ng th e ACSIS pipel i ne (ICavanagh et al.l l2008l ) in a similar 


manner as described in Currie (2013) and Sadavov et al. ( 2013l i. In brief, the HARP data were 
reduced using a two-stage, fully automated methodology. First, we conducted a quality assurance, 
which mostly removed those spectra degraded by a variety of instrumental effects. Second, we 
performed an iterative procedure that used spectral cubes to estimate linear baselines and to 
identify emission features to be masked for a more refined baseline subtraction. For the Bl-F data, 
one iteration was sufficient to characterize the baselines. Detailed descriptions of the methods used 
can be found in Jenness et al. (2015, submitted). 


Finally, we used a main-beam efficiency of rjMB = 0.61 for our reduced spectra. The effective 
beam of these data is 17.3" , which we convolved with a 31.9" beam to get our data at an equivalent 
resolution of 36.3", or the resolution of our 500 /im Herschel observations. The 1 a line sensitivities 
at this spatial resolution are ~ 0.2 K and ~ 0.3 K for the ^^CO (3-2) and C^®0 (3-2) data, 
respectively. 


2.3. Submillimeter Telescope 

The Bl-F clump was mapped using the 10 m of the Arizona Radio Observatory (ARO) Sub¬ 
millimeter Telescope (SMT) in dual polarization, 4-IF mode over ~ 35 hours in March 2012. We 
used the 1 mm ALMA prototype receivers to make position-switching on-the-fly (OTF) maps with 
the upper sideband tuned to (2-1) at 230.538 GHz and the lower sideband tuned to 

(2-1) at 220.3987 GHz with a ~ 5.0 GHz IF. With the 250 kHz filterbank backends, the velocity 
resolution is ~ 0.34 km s“^ and the spatial resolution is ~ 33.7" at 220 GHz. At this time, we only 
discuss the ^^GO (2-1) data. 

Point and focus observations were made on Venus or Jupiter. Line calibrations to measure 
the sideband spillover were made using the standard SMT line calibrator in Orion A (J2000.0: 
05:35:14.48, -05:22:27.6). The beam efficiency (ry « 0.62) was measured using the average efficiency 
from observations of Venus and Jupiter. The beam efficiency was stable (within 5%) throughout 
the observations. 

To avoid artifacts introduced by long scans, we divided Bl-E into smaller subregions, mapping 
each block separately. We used four overlapping subregions, where each subregion was 9' x 7', 
with a 2' overlap, for a total map area of ~ 190 arcmin^. We initially used an OTF scanning rate 
of 10"/s, but switched to 20"/s to speed up the mapping. Each region was observed in orthogonal 
scans (along RA or deck) with 10" spacings at least once at 10"/s or twice at 20"/s along each 
direction. We selected the coordinates 03:35:31.0, 31:45:31.0 (J2000.0) for our off-position. 
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The observations were reduced using the gildas reduction software, CLASS. We reduced 
the horizontal and vertical polarizations separately using an initial single-order baseline fit. The 
telescope was fairly stable over the observations, and we rejected only two scans in one subregion 
due to unusual spectral signatures. The spectral line data were combined onto a evenly spaced grid 
and smoothed to 36.3" resolution, the resolution at the 500 /rm beam. Several rows of spectra had 
curved baselines. To improve these baselines, we subtracted an additional quadratic function using 
a wide velocity range (-11 km s“^ to 28 km s“^). The final 1 a rms of the smoothed maps is ~ 0.05 
K towards the center (highest sampling) and ~ 0.15 K towards the corners of the map. 


2.4. IRAM 30 m Telescope 

We also fully mapped Bl-E with the IRAM 30 m telescope in August 2012 using the EMIR 
receiver and VESPA backend to observe simultaneously (1-0) at 110.2014 GHz, (1-0) at 
109.7822 GHz, and (2-1) at 219.5604 GHz. We made frequency-switching OTF maps with a 

throw of 3.9 MHz in the 90 GHz band and 11.7 MHz in the 230 GHz band. The respective velocity 
resolutions are ~ 0.106 km s“^ and ~ 0.053 km s“^ and the respective spatial resolutions are ~ 24" 
and ~ 12" at 110 GHz and 220 GHz. 

Like our SMT observations, we mapped Bl-E with orthogonal scans in four smaller, overlapping 
subregions for a total coverage of ~ 190 arcmin^ in all three tracers (see Section 12.31) with a step 
size between scans of 5" and a scan rate of 7.5"/sec. Point and focus observations were done using 
Uranus after full scans of each subregion. Line calibrators were observed while Bl-E was in transit. 
We adopted beam efficiencies and forward efficiencies of r/fegam = 0.78 and Eg// = 0.94 at 110 GHz 
and ribeam = 0.61 and F^ff = 0.93 at 220 GHz. The total on sky time was ~ 13 h. 

The IRAM 30 m observations were reduced using CLASS. We removed an initial single-order 
baseline from all spectra before folding the spectra. For the G^®0 (1-0) and G^®0 (2-1) observations, 
we subtracted additional single-order or second-order baselines to remove large-scale ripples after 
the spectra were folded. As with the previous line data, we smoothed the data to a similar grid as 
our the 500 fim Herschel data. The hnal 1 a rms of the smoothed maps range from ~ 0.03 — 0.1 
K for (1-0), ~ 0.03 - 0.1 K for G^^O (1-0), and ~ 0.1 - 0.2 K for G^^O (2-1) between the 

center of the map (highest sampling) and the corners (lowest sampling). 


3. Results 


Figured] shows the Herschel-deiived column density map for the Bl-E clump. Golumn densities 
(and line-of-sight dust temper atures) were de t ermin ed from modihed blackbody fits to the Herschel 
data in the same manner as ISadavov et al.l (|2014l ). where we assumed a fixed dust opacity law 


corresponding to /3 = 2 and kq = 0.1 cm^ g ^ at 300 /rm, which is the adopted dust opacity law by 


the HGBS consortium (e.g., see lRov et al.ll2014l ). and we applied average colour corrections to the 
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Herschel field. Moreove r, we assume d a mean mole c ular w eight per hydrogen molecule of ^ = 2.8 
([Kauffmann et al.l 120081 1. Note that ISadavov et al.l (1201211 adopted a mean m olecular weight per 
unit mass (/i = 2.33). Additionally, the observations in ISadavov et al.l (|2012l l were reduced using 
HIPE 5, and due to updated calibration references, the PACS fluxes changed by < 20% between 
HIPE 5 and HIPE 6 (B. Altieri 2013, private com munication). Thus, th e column densities presented 
here are lower by ~ 32% compared to those from ISadavov et al.l (j2012l ). 
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Eig. 1.— Herschel-derived column density map for Bl-E based on modified black b ody fi ts to the 
160 — 500 ^m SEDs. The position of the nine substructures from ISadavov et al.l (120121 ) are also 
given, with circles corresponding to the map resolution (36.3''). The grey contours correspond to 


l21 2 


column density levels of 7.0 x 10^^ cm" ^,8.5x10^^ cm and 10.5x10 
contour represents the outer boundary of the clump (N(H2) = 5.0 x 10^^ cm 
small crosses correspond to spectra shown in Figure [2l 


and the thick black 
Av ~ 7 mag). The 


Figure m also gives the positions of the nine Bl-E substructures identified bv ISadavov et al 


(j2012l . see paper for coordinates). These substructures are typically < 0.1 pc in size with masses 
of < 1 Mq, and a verage densitie s of ^ 10^ cm“^. For more physical information about the Bl-E 
substructures, see ISadavov et al.l (120121 ). 


Table [T] lists the mass, approximate angular dimensions, area, effective radius, and median 


















































- 7 - 


temperature of the Bl-E clump using the Herschel data alone. These values correspond only to the 
region of Bl-E with N(H2) > 5 x 10^^ cm“^ (thick black contour in Eigure[T]). Mass was measured 
from the column densities, clump area was determined from the number of pixels above the column 
density threshold, and the effec tive radius was me asured assuming Bl-E is circular, Reff = \/ ^/tt, 
assuming a distance of 235 pc (jHirota et al.ll2008l b The median temperature comes from the line- 
of-sight dust temperatures determined from the modihed blackbody fits to the Herschel data. 


Table 1: Herschel-deiived Bl-E Properties 


M (Mo) 

Angular Size 

A (pc^) 

Reff (pc) T(K) 

88 

15.5' X 10.5' 

0.57 

0.43 13.8 


3.1. Velocity Structure 


Eigure [2] shows three sample IRAM (1-0) and (1-0) spectra for the sample pixels 

toward the north-west (NW) corner, center (C), and south-east (SE) corner of the Bl-E clump (see 
Eigure d]). The IRAM (1-0) line emission has a main velocity component at ~ 7.5 km s“^ 

with less prominent emission at ~ 3 km s“^. The integrated intensity associated with the ~ 7.5 
km s“^ component dominates over the ~ 3 km s“^ component by a factor of ~ 5 on average. The 
C^®0 (1-0) observations barely show any emission at ~ 3 km s“^, indicating that the gas associated 
with the 3 km s~ ^ component is more diffuse than the gas found at the ~ 7.5 km s“^ (see also, 
Sadavov et al.ll2012l b Hereafter, we will only discuss the 7.5 km s“^ velocity component, unless 


stated otherwise. 


The IRAM (1-0) and IRAM C^®0 (1-0) spectra at ~ 7.5 km s“^ reveal very complicated, 
non-Gaussian profiles indicative of several sub-components at very similar, overlapping velocities. 
The best-fit velocity profiles in Figure [2] show single Gaussian fits for each of the 3 km s“^ and 7.5 
km s“^ components, whereas many IRAM (1-0) spectra require five unique velocity subcom¬ 

ponents to fit well the observed velocity profiles at ~ 7.5 km s“^, and similarly, the IRAM G^®0 
(1-0) spectra need three components at ~ 7.5 km s“^. The SMT ^^GO (2-1), JGMT ^^GO (3-2), 
and IRAM G^^O (2-1) line emission also show complex line prohles. 

Figures [3] and m show respectively the integrated intensity maps for all the ^^GO and G^®0 
transitions across the main velocity structure over 5—9 km s“^. Both figures include N(H2) contours 
to highlight the regions with higher column densities of dust. We find that the ^^GO observations 
show some agreement with the dust, whereas the IRAM G^®0 (1-0) and (2-1) integrated intensities 
show excellent correspondence with the dust-derived column densities. Although the right panel of 
Figure 0] hints at possible JGMT G^^O (3-2) emission, particularly towards the Bl-El region (see 
Figured]), any detections are very weak (i.e., < 3 a). We do not use the JGMT G^®0 (3-2) data in 
our subsequent analyses. 
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Fig. 2.— Sample IRAM spectra for (left) (1-0) and (right) C^®0 (1-0) towards the north-west 
corner (NW), center (Cen), and south-east corner (SE) of Bl-E. The pixels used for the spectra 
are shown in EigurelH The red curves represent the best-fit line profiles using a single Gaussian fit 
for the 3 km s“^ component and 7.5 km s“^ component, respectively. The dotted grey lines show 
the 3 a rms level. 


3.2. Velocity Gradients 


The IRAM (1-0) and IRAM (1-0) spectra show a radial velocity gradient, where the 
centroid velocities increase from the NW to the SE (see Figure [2]). Figure [5] shows position-velocity 
diagrams for the IRAM (1-0) and IRAM (1-0) line transitions. We find consistent 

flux-weighted gradients for both line transitions of 0.94 km s“^ pc“^ from the IRAM (1-0) 

observations and 1.0 km s“^ pc“^ for the IRAM (1-0) observations at an a ngle of ^ 20 — 30° 


north-to-west. This gradient of ~ 1 km s ^ pc ^ is consistent with the results of iKirk et al.l (j2010l ) 


for the Bl-E clump. We also find a similar gradient and angle with the 3 km s ^ component. 

On much larger scal es, a steep velocity gradient across the Perseus cloud has been well identified 
in previous studies (e.g.. lRidge et al.ll2006l : ISun et al.ll2006l) . Figure [6] compares the channel maps of 
Bl-E to those of L145 5 (clump west of B l-E) and IC348 (clump east of Bl-E) using data from the 
COMPLETE survey (jRidge et al.ll2006l ). The line-of-sight velocities in Perseus typically increase 
from the south-west (SW) to the north-east (NE), as illustrated by the channel maps of L1455 and 
IC348, as well as their centroid velocities. Moreover, the COMPLETE ^^CO (1-0) and ^^CO (1-0) 
data give very similar gradients. 
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a (2000) 


Fig. 3.— Integrated intensity maps for line emission over the velocity range of 5 — 9 km s“^ 

only. The left panel shows IRAM (1-0), the middle panel shows SMT (2-1), and the 

right panel shows JCMT (3-2). The lower greyscale limit corresponds to the 3 a uncertainty 

in the integrated intensity. The contours show Rersc/ie/-derived column density levels of 5.0 x 

and 10.5 x 10^^ cm“^ from the RGBS data (see 


10^^ cm 7.0 X 10 ^^ cm ^,8.5 x 10^^ cm ^ 


Sadavov et al.ll2014l l. 
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Fig. 4.— Same as Figure [3l but showing the IRAM (1-0), IRAM (2-1), and JCMT 

(3-2) integrated intensities. 


To characterize better the velocity gra dients elsewhere i n Perseus, we used the (1-0) 

observations from the COMPLETE survey (jRidge et al.ll2006l i in a similar manner as our IRAM 
observations of Bl-E. We broke up the Perseus cloud in to rectangular region s identified by-eye to 
enclose B5, IC348, Bl-E, Bl, L1455, and L1448 (e.g., see Sadavov et al. 2014), as well as addi tional 
regions encompassing the “bridge” material between IC348 and Bl-E (see lBallv et al.ll2008l ). We 
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Fig. 5.— Position-velocity diagram for (top) IRAM (1-0) and (bottom) IRAM C^®0 (1-0). 

Lines were averaged across strips in declination (left) or Right Ascension (right). For (1-0), 

contour levels correspond to average Tmb values of 0.5 K, 1.0 K, 1.5 K, 2.0 K, 4.0 K, 5.0 K, and 
5.7 K. For (1-0), contours correspond to average Tmb values of 0.04 K, 0.08 K, 0.2 K, 0.5 

K, and 1.0 K. The red lines correspond to the weighted mean velocity for Tmb > 1-0 K (^^CO) or 
Tmb > 0.1 K (Ci^O). 


excluded NGC1333 due to its large sample of outflows (IHatchell et al.ll2007l : lArce et al.ll2010lL Note 
that we generally selected more material than iKirk et ah ( 2010 ) in measuring the gradients. Figure 
[7] shows the radial velocity gradient orientations across Perseus using the same technique outlined 
in Figure [5l Most regions have gradients of ~ 0.3 — 0.6 km s“^ pc~^, where the centroid velocities 
increase from the SW to the NE at position angles of ~ 30 — 50° north-to-east. Thus, the velocity 
gradients toward Bl-E and the nearby “bridge” region stand out significantly for increasing instead 
in directions perpendicular to the velocity gradients found in the other Perseus clumps. 
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Fig. 6.— Channel maps of (1-0) observations of IC348 (left), Bl-E (i niddle), and L145 5 

(right). The IC348 and L1455 maps use data from the COMPLETE survey ( Ridge et al. 2006 1. 
whereas the Bl-E data use the IRAM ^^CO (1-0) data presented here. The contours for IC348 and 
L1455 correspond to 20%, 40%, 60%, and 80% of the peak line brightness of each channel. The 
contours for Bl-E correspond to 50-90% of the peak line brightness (in steps of 10%). The observed 
velocity channel is given above each map. 


3.3. Turbulence 

To quantify the level of turbulent motion in Bl-E, we measured the non-thermal velocity 
dispersion for both IRAM ^^CO (1-0) and IRAM C^®0 (1-0) line emission using two methods; (1) 
we adopted the median velocity dispersion from the single-Gaussian fits (see Section [3T]) : (2) we 
summed all the spectra, where each spectrum was re-centered prior to stacking them (to correct 
for the observed velocity gradient; see previous section). Table [2] gives the velocity dispersions 
associated with each method and each molecule. For both molecular tracers, the velocity dispersion 
is slightly higher when using the stacked spectra compared to the median velocity dispersion of all 
spectra, though both methods agree within uncertainties. Additionally, our line widths for Bl-E 
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Fig. 7.— Velocity g radient orientation s across Perseus as determined from (1-0) observations 

from COMPLETE (jBidge et al.l 120061 ). Background image shows the integrated intensity map for 
(1-0) using the COMPLETE data. Arrows show the direction of centroid velocities from bluer 
values to redder values (they are not scaled to the gradient magnitude). We excluded NGC1333 
from this analysis (see text). 


agree well with those found in other (more active) clumps of Perseus (e.g., iHatchell et al.ll2005l : 
Kirk et al. 2010). 


We estimated the thermal and non-thermal line contributions using the median dust temper¬ 
ature of ~ 14 K, as determined from modified blackbody fitting to the Herschel SEDs (see also 
Appendix]^, as a proxy for the gas kinetic temperature. Assuming an equivalent gas kinetic 
temperature of 14 K, the thermal line widths of ^^CO and C^®0 are ~ 0.06 km s“^, whereas the 
observed velocity dispersions are a > 0.3 km s“^. Although the gas and dust are likely decoupled 
at the densities th e ^^CO and C^^O l ine emission trace, th ese temperatures should agree within a 
factor of two le.g.. lYoung et al.ll2004l : ICeccarelli et al.ll2007l L Thus, the expected thermal velocities 
will be < 0.1 km s“^ (e.g., for Tk < 2Tdust) and the ^^CO and C^®0 spectra cannot be explained 
by thermal broadening alone. Rather, we assume the observed line widths represent the turbulent 
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Table 2: Turbulent Velocity Dispersions Across Bl-E 


Line 

Median'^ 

Stacked Spectrum*^ 

i^GO ( 1 - 0 ) 

0.7 km s“^ 

0.8 km s“^ 

0^*0 ( 1 - 0 ) 

0.4 km s“^ 

0.5 km s“^ 


“The distribution of velocity dispersions across Bl-E are non-Gaussian. For the (1-0) and C^®0 (1-0) data at 

N(H 2 ) > 5.0 X 10^^ cm“^, ~ 75 % of spectra have velocity dispersions within 0.13 km s“^ and 0.15 km s“^ of the 
respective median velocity dispersion. 

^’Single Gaussian fits yield robust velocity dispersions with uncertainties better than 1%. 

(non-thermal) velocity dispersion. 


3.4. Substructures 


Sadavov et all (|2012l ) identified nine substructures in Bl-E (see also, Figured]) and character¬ 


ized them with NH 3 (1,1) emission. Figures [ 8 | and [9| compare the (1-0) and (2-1) spectra of both 
and C^®0, respectively, with the NH 3 (1,1) centroid velocities (dashed line) of these nine 
objects. The NH 3 (1,1) data trace the denser central material of each substructure, whereas the 
and data trace better the source envelopes and more diffuse clump, respectively. In 

general, most substructures show a good agreement between the N H 3 (1,1) centroid ve l ocities and 
the velocity peak s , in agreement w ith previous studies (e.g., IWalsh et al.l l2004l : iKirk et al 


20071 : 1 Andre et al.l 120071 : iKirk et al.l 120101 ) . This similarity supports the notion that th e dense ma¬ 
terial within cores forms quiescently within their envelopes (see also. lOffner et al.ll2008l ). The 
spectra show more deviations in centroid velocity (and higher line widths), which is unsurprising 
as is a better tracer of lower density (~ 10 ^ cm“^) gas. 


The Bl-El and B1-E4 substructures are the only objects with significant (> 3 a) deviations 
between the peak-line velocities of and NH 3 (Figures [ 8 | and [9|). In particular, these objects 

have asymmetric profiles with a narrow Gaussian-like peak and a broad blue-shifted “knee,” 

where the NH 3 emission peaks between the Gaussian-peak and the “knee.” In addition to these, 
the B1-E6 source has a double peaked G^®0 spectrum (a feature seen in both transitions) with 
single-Gaussian NH 3 (1,1) emission at velocities directly between the two peaks. Normally, a 

double-peaked spectrum with brighter emission toward the bluer side would suggest self-absorption 
and perhaps infall. Nevertheless, emission is generally optically thin (see Section Hj) and as 

such, will not produce such a feature. Moreover, the spectra, which should be more optically 

thick, do not show this double-peaked profile. Thus, the B1-E6 spectra may correspond to multiple 
velocity subcomponents along the line of sight. For example, the bluer velocity peak matches 
the centroid velocity of the nearby B1-E7 source, whereas the redder peak appears associated with 
a separate velocity structure east of the B1-E6 and B1-E7 ridge. We discuss the possibility of 
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Fig. 8 .— Comparison of the IRAM (1-0) and IRAM (1-0) line profiles for the nine 

densest Bl-E substruct ures. The dotted vertical line corresponds to the respective NH 3 (1,1) 


centroid velocity from ISadavov et al.l (j2012l i . 


multiple subcomponents in Section [5l 

In addition to similar core-envelope velocities, several studies have compared the core-to- 
core centroid velocit y dispersion with the turbulence in their surrounding environment. From 
Sadavov et al.l (j2012l ) , the substructures have an NHa-derived core-to-core relative motion of cTsys ~ 
0.24 km s“^. This core-to-core centroid velocity dispersion matches well the dispersion of the 
(1-0) peak velocities (~ 0.29 km s“^), but is factors o f 2-3 lower t han t he typical velocity disper¬ 
sions observed across the entire clump (see Table m- iKirk et al.l ((20101) similarly found that the 
starless (and protostellar) cores in other regions of Perseus systematically had much lower rela- 
tive motions tha n the turbulence detected in the more diffuse, ambient clump material (see also, 
Kirk et al.ll2009l i. Therefore, these core precursors do not appear to move in a manner similar to 
the Bl-E clump as a whole, but rather, they have more homogenous velocities like those seen in 
more evolved dense cores. We also find that the NH 3 detections show no indication of a velocity 
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Fig. 9.— Same as Figure [ 8 ] but for the SMT (2-1) and IRAM (2-1) transitions. 


gradient, whereas the C ^°0 data clearly show a ~ 1 km s“^ pc“^ radial velocity gradient across 
Bl-E. Thus, the Bl-E substructures appear to be decoupled (i.e., following a distinctly different 
velocity distribution) from the moderately dense gas as traced by the emission, suggesting 

that the systemic motions of dense cores somehow homogenize very quickly in their evolution (see 
also. Section [5.2.2p . 

The Bl-E substructures have very broa d NH.-^ emission c o mpare d to typical observations of 
dense gas in prestellar cores. For example, iRosolowskv et al.l ([2003) conducted a targeted NH 3 
survey of Perseus, and from their results, the prestellar cores have NH 3 velocity dispersions a < 0.25 
km s“^, with an median value of a ~ 0.15 km s“^. In comparison ^ most of the Bl-E s ubstructures 
have much broader NH 3 velocity dispersions of cr > 0.3 km s“^ ( Sadavov et al.ll2012l ). The most 
line-broadened substructures have similar velocity dispersions (~ 0.4 km s“^) as the larger-scale 
gas. If these broad line widths represent the turbulent motions within these objects, the 
Bl-E substruc tures have significant internal pressures. Assuming a typical substructure density of 


10^ cm ^ (see ISadavov et al.ll2012h . the substructures have internal pressures of Psub/k = pa^ = 
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(0.3 — 6) X 10^ K cm In comparison, the volume-a veraged gravitational pressure of the Bl-E 
clump can be estimated from iBertoldi &: McKeel (|l992l l as 


P f M 


R 

1 pc 


-4 


K cm 


-3 


( 1 ) 


where M is the clump mass, R is the observed clump size, ai is a scaling factor corresponding to a 
nonuniform density distribution, and ((t> n) is a scaling fact or corresponding to the cloud geometry 
{4>g = 1 for a perfect sphere). Following iLada et al.l (|2008l l. we assume ai = 1.3, appropriate for a 
self -gravitating cloud. We a lso adopt {4>g) ~ 1.01, since Bl-E is fairly circular (~ 1.2 aspect ratio; 


see iBertoldi &: McKeel Il992l ). Thus, the Bl-E volume-average gravitational pressure is ~ 3 X 10^ 
K cm“^, which is very similar to the internal pressures for the Bl-E substructures. Moreover, if 
the substructures are centr ally concentrated, their surface pressures may be much lower than their 
average internal pressures (jLada et al.l l2008l l . suggesting that the external pressure of the clump 
is sufficient to keep the Bl-E substructures confined. Such pressu re-confined stru ctures have been 
seen more recently in different environments, including Taurus (jSeo et al.l l2015l l and Ophiuchus 
( Battle et al. 2015). 


4. Radiative Transfer Models 


We used the radiative transfer code, RADEX, by Ivan der Tak et al.l (120071 ) to characterize 
the abundances of and C^®0. In brief, RADEX is a non-LTE radiative transfer code that 

models molecular transitions to determine the expected line strengths, optical depths, and energy 
level populations for different transitions of a given molecule. The code assumes a homogenous, 
spherical cloud with a user-determined temperature and abundance. Since our spectra contain 
multiple velocity components, we used the results from our single-Gaussian fits (see Figure [2]) to 
determine the observed line brightnesses and widths for all our transitions. Thus, our RADEX 
results represent the bulk properties of the Bl-E clump. 

We ran RADEX for individual pixels across our and C^®0 observations. (Recall that all 

data have been convolved to a common resolution and rebinned to a common grid.) We selected 
only those pixels with N(H 2 ) > 5 x 10^^ cm“^, to ensure that we had strong detections > 5 cr for 
all transitions considered. For each pixel, we adopted a background temperature of 2.73 K and the 
line widths from our single Gaussian fits. Additionally, we assumed the kinetic gas temperature is 


equa. 1 to the observed dust temperature from SED fitting of the Herschel data (see ISadavov et al 


20141 ). Since the Bl-E is a moderately dense region, the gas and dust temperatures may not be 


coupled. We expect differences of only a few Kelvin, which produce relatively minor uncertainties 
in the models (see Appendix lAl) . 


We applied 2500 combinations for the H 2 number density, n(H 2 ), and the molecular column 
density, N(^^GO) or N(C^®0), in our analysis. For n(H 2 ), we adopted values between 10^ cm“^ < 
ii(H 2 ) < 10^ cm“^ in equal spacings of Alog[n(H 2 )] = 0.04. For the molecular column densities, 
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we adopted 5 x cm“^ < N(^^CO) < 5 x 10^® cm“^ and cm“^ < N(C^®0) < 10^® cm“^ 
with equal spacings of Alog[N(X)] = 0.04. For each density/column density pair, we compared 
the model line brightnesses and the observed line brightnesses for all observed transitions, and we 
adopted the parameters associated with the best-fit model. 

Figure fTOl shows the best-fit column densities for and from the RADEX models as 

well as their corresponding X-factors, using our Herschel-dex'wed H 2 column densities (see Figure 
d]) to determine N(X)/N(H 2 ) for both and C^^O. We find that the column densities 

trace the densest material (see contours) very well, whereas the column densities do not follow 
the dense material as closely. For the emission, we find moderate optical depths of ~ 3 at 

all J-lines connecting the material around B1-E2 to B1-E8 and B1-E9. On average, the J = 2 — 1 
transition has the highest optical depths, whereas the J = 1 — 0 has the lowest optical depths. In 
contrast, the emission is optically thin, suggesting that we are tracing all the emission along 

the line of sight. 


Eigure 1111 compares the observed valu es of N(C^^O ) and N(H 2 ) with the X(C^^O) relations 
from Frerking et al. (1982, for Ophiuchusfl IPuvert et ^ ( 1986 . for Taurus), and Kainulainen et al 


(200g, for Chameleon). The Frerking et al. relation generally underestimates the abundances 
and the Duvert et al. relation generally overestimates them. The Kainulainen et al. relation agrees 
well with N(C^®0) at low column densities (e.g., N(H 2 ) < 8 x 10^^ cm“^), but underestimates 
the abundances at higher column densities. Moreover, for N(H 2 ) < 8.0 x 10^^ cm“^, we 

find that the column densities are much lower, on average than what is expected by the 

X(C^®0) relations. These lower abundances could be due to selective dissociation of C^®0 

molecules, which de pends greatly on the se l f-shiel ding of the molecule (abundance) and geometry of 
the clump. Indeed, Ivan Dishoeck &: BlackI ( 19881 ) showed that the photodissociation rate of C^®0 


in the centers of clouds is generally much higher than the photodissociation rate of and 
due to self-shielding. At the very highest column densities (N(H 2 ) > 1.0 x 10^^ cm“^), there is far 
less scatter in the abundances, suggesting that the gas is more shielded. 


To account for selective photodissociation better, we used a linear least-squares regression 
to fit simultaneously two distinct linear relations to the lower-density material (N(H 2 ) ^ 1 x 
10^^ cm“^) and to the higher density material (N(H 2 ) ^ 1 x 10^^ cm“^), excluding the data 
points associated with B1-E2 (black symbols). We find a best-fit “break” in the two relations (e.g.. 


where becomes well shielded) at N(H 2 ) = (10.7 ± 0.2) x 10 ^^ cm ^ {Ay 15 mag), with 

X(C^®0) ~ 2.7 X 10“^ at column densi ties below this break valu e and X(C^^O) 2 .2 x 1 0~^ at 

higher column densities. Similarly, both iFrerking et al.l (jl982l ) and iKainulainen et al.l (j2006|) found 
that the N(C^®0) relation tends to flatten at higher Ay values, though the break identified by 
Frerking et al. occurs at a column density below what is considered here for Bl-E, e.g., at Ay ~ 4 
mag, whereas we only use data at an equivalent of Ay > 7 mag. Therefore, self-shielding 

may start at Ay ~ 4 mag, resulting in a first break in X(C^®0), and then, depending on the cloud 


^We use only the high-Av relation of Frerking et al. to better match our Bl-E data with Ay > 7 mag. 
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Fig. 10.— Top: Column density maps of (left) and (right) from our best-fit RADEX 

models. Bottom: Ratios of (left) and (right) column densities with the Herschel-derived 

H 2 column density from the RGBS data. Contours correspond to N(H 2 ) = 7.0 x 10^^ cm“^,8.5 x 
10^^ cm“^, and 10.5 x 10^^ cm“^. 


geometry, the C^®0 gas is mostly self-shielded by Ay ~ 15 mag. We note that this higher Ay 
break for well-shielded gas is similar to the Ay ~ 12 mag break iKainulainen et ahl (j200m found for 
Chameleon. 


Our best-ft relation for X(C^^O) excluded the column densities associated with B1-E2, which 
appear systematically lower than all the other Bl-E substructures (see Eigure [TT]l . This deviation 
could indicate that the B1-E2 source has different dust properties than observed elsewhere in Bl- 
E. Eor example, if the true dust opacity toward B1-E2 is higher by a factor of two (e.g., due to 
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Fig. 11.— Comparison between N(C^^O) and N(H 2 ) for Bl-E. The black points represent the 
column densities assoc iated with B1-E2. Th e red curve shows the X(C^®0) relation for Ophiuchus 


from 

Duvert et al. 

(198611 

Kainulainen et al. 

(20061 


ii 8 r 


Since all three studies primarily fit data w ith lower Ay values (e.g., 


Ay < 7 mag), we used the N(H 2 )/^y relation from Bohlin et al. ( IQTsl i of N(H 2 )/^y = 9.4 x 10^° 


cm ^ mag The Bohlin et al. ratio better rep resents l ower c olumn density environments (e.g., 

^20 


cm than the relation from IPraind ( 20031 1 of N(H 2 )/Ay = 6.9 x 10 


N(H2) < 6 X 10^ 

mag“^ ( Rov et al.ll2014l ). The green curve corresponds to our best-fit relation (excluding the B1-E2 


data points) from a two-component linear least-squares regression. We find a best-fit break in the 
curve at N(H 2 ) = (10.7 ±0.2) x 10^^ cm 

than the break and X(C^®0) ~ 2.2 X 10 ^ at higher column densities. 


,21 2 ^ with X(C^®0) ~ 2.7 X 10 ^ at lower column densities 


dust coagulation, see lOssenkopf fc Hennindll994l b then its X(C^®0) ratios would match the rest of 
the substructures. Nevertheless, to produce such a significant difference in the local dust opacity 
toward just B1-E2, the physical p roperties of this object must be significantly different than the 
other substructures. Eor example, lOssenkopf fc Hennind (119941 1 find that densities that varied by 
two orders of magnitude will have dust opacities (from coagulation alone) that differ by less than 
a factor of two in the Herschel ban ds. Since B1-E2 has a similar density and dust temperature as 
the rest of the substructures (e.g., ISadavov et al.ll2012l: Pezzuto et al. in preparation), we would 
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not expect a correspondingly large increase in dust opacity. 

In addition to the lower abundances, B1-E2 has an enriched chemistry of ni trogen-bearing 


mole cules and carbon-chain molecules compared to the other substructures (see ISadavov et al 


20121 ). indicative of a distinct chemistry toward only this substructure. One explanation for this 


distinct chemistry is that the carbon-bearing molecules are more depleted from the gas phase for Bl- 
E2 and that the lower X(C^®0) factors for B1-E2 are due to this freeze-out of C^®0 molecules. Many 


Bacmann et al. 

2009: 

Duarte-Cabral et al. 

2010: 

Hernandez et al. 

2011 

) and 

Bergin et al. 


(|2002l l found a flattened, slightly decreasing profile in both C^°0 and C^'O emission toward very 


high column densities in B68, in good agreement with our flattened N(C^^O) profile for B1-E2. We 
discuss depletion from molecular freeze-out in more detail in Section [5.2.11 


5. Discussion 


5.1. Formation of High Density Objects 

The formation mechanism for the Bl-E substructures must explain why the Bl-E clump is 
only condensing to form cores now, whereas the other clumps in P erseus have formed populations 


of dense cores and YSOs (e.g., iKirk et al.M2006l : lEvans et al.l 120091 '). One simple explanation could 


be that the condensation of the Bl-E clump itself was triggered only recently. Bl-E has a radial 
velocity gradient nearly perpendicular (in the plane of the sky) to the typical velocity gradient 
associated with the main Perseus cloud (see Section . Thus, this clump may be impacted by 
large-scale processes such as shearing flows or collisions that have only occurred recently. Indeed, 
we see unusual, non-Gaussian spectra t oward Bl-El and B1-E4 t hat match well the spectra 


associated with colliding clouds modeled by iKeto fc Lattanzid (119891 ). Since these two sources 


are located along the western-edge of Bl-E at the junction between the Bl-E velocity gradient 
and the Perseus overall gradient, they are most likely to be affected by the interaction of these 
gas flows (see Section [3lT)l . Nevertheless, li 
motions and filamentary mo r phologies (e.g . 


20041 : IClark &: Bonnelll l2005l : iKlessen et al.l [ 2003 ; 


Klessen et al. 

2000: 

Bate et al. 

o 

o 

Li &: Nakamura 


Basu et al. 


2009aj), and the Bl-E clump and its 


substructures shows no evidence of these properties. 

The relationship between the Bl-E velocity gradient and that of the Perseus cloud remains 
unclear. For example, perpendicular velocity gradients may arise from (1) independent clouds along 
the line of s ight, (2) material flowing into or out of the system, (3) dynamic turbulence on small 
scales le.g.. [Smith et al.ll2012h . or (4) strong magnetic fields. At this time, we cannot distinguish 
between these possibilities to explain the perpendicular velocity gradient between Bl-E and the 
Perseus cloud. Specific predictive simulations of clump fragmentation triggered by low-velocity 
shearing (or colliding) flows are needed to determine the origin and impact of the perpendicular 
velocity gradient in Bl-E. 
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If the condensation of the Bl-E clump was not triggered by an external proce ss, its star 
formation could be delayed compared to the other Perseus clumps. ISadavov et al.l (j 2012 l ') suggested 
that Bl-E may be under the i nfluence of a strong, loc alized magnetic field. Bl-E has uniquely strong 
R-band polarization vectors ([Goodman et al.lll990l l. and a strong magnetic field has been shown 
to both delay star formation and to form objects in a l oose configuration without filaments as seen 
in Bl-E (e.g.. lLi fc Nakamural 12004 : iBasu et al.ll2nn9al ). Alternatively, Bl-E and its su bstructures 


may b e influenced by significant turbulent processes (see Sections 13.11 and l3.4p . Indeed. IPon et al 
( 2014l l found excess emission in high-J lines toward the B1-E5 substructure, which they 

interpreted as evidence of low-velocity shocked gas. Such small-scale turbulent shocks are expected 
to form few so urces in isolated mo rphologies on longer time scales than large-scale, high velocity 
processes (e.g., iKlessen et al.l 120001 1 . Thus, the loose configuration of the Bl-E substructures and 
their delayed formation with respect to the other Perseus clumps may indicate that this region is 
dominated by small-scale turbulent shocks or magnetic fields. 


5.2. Evolution of High Density Objects 

5.2.1. Comparisons with More Evolved Perseus Clumps 

The Bl-E clump and its substructures have some similarities and differences compared to 
other, more evolved regions of Perseus. For example, we find typical line widths of a ~ 0.8 
km s“^ and a ~ 0.4 km s“^ from the IRAM (1-0) and IRAM (1-0) observations, 

respectively, in excellent a g reeme nt wi th the equivalent tur bulent results for the Perseus cloud 
in general from iKirk et al.l (|2010l l and iHatchell et al.l ()2005l l. Moreover, the Bl-E substructures 
appear decoupled from their parent clump. They do not follow the velocity gradient seen in the 
larger -scale and emission and their central densest material, as traced by NH 3 emission 

(from ISadavov et al.ll 2012 l l. have more hon iogenous ve l ocitie s compared to the moderately dense gas 
traced by either or C^®0. Similarly. iKirk et al.l ( 2010l ) found that the dense (prestellar) cores 

in the more evolved regions of Perseus also had more homogenous centroid velocities compared to 
the surrounding, larger-scale gas. Therefore, Bl-E has global kinematic properties similar to those 
of the other Perseus clumps. 

In contrast, the Bl-E substructures were not well detected by SCUBA or Bolocam, whereas 
both instruments were able to detect well prestellar cores in other regions of Perseus. Additionally, 
the Bl-E objects have broad, supersonic NH 3 line widths (a > 0.3 km s ~^), whereas the genera l 


prestellar core population in Perseus are more quiescent (a < 0.25 km s ^: IR osolowskv et al.ll2008l ) 


Such broad line widths could indicate either significant in ternal turbulent motions or multiple , 
unresolved subcomponents within the GBT beam (e.g., see lSmith et al.ll2012l : iBailev et al.ll2015l ). 
Since the Bl-E clu mp is not a very dense region and NH 3 (1,1) is not well detected off the column 
density peaks Isee iRosolowskv et al.ll2008l l. line-of-sight coincidences of unrelated subcomponents 
of dense gas are unlikely to explain the broad NH 3 (1,1) emission for all the substructures (except 
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possibly B1-E6, where we detect two distinct velocity components in the line emission). 

Instead, the central dense gas itself may be fragmented into different velocity subcomponents. 
Without higher resolution observations, however, we cannot identify multiple sources. 

One of the Bl-E substructures stands out as unique. B1-E2 is the only object with transonic 
NH 3 (1,1) line widths (u ~ 0.1 km s“^), similar to what is detected in typical prestellar cores. 
Moreover, B1-E2 also has relatively lower column densities than the other substructures for 

the same column densities (see Figure n]), which could indicate depletion uniquely toward 

this object, another property associated with prestellar cores. This decrease in abundance 

is well within the uncertainties of X(C^®0), however. F or example, N(H^) is uncert ain within a 
factor of 2 from assuming a single dust opacity (e.g., see lOssenkopf &: Hennindll994l ). Moreover, 
our RADEX models assume a single system temperature, single system velocity, and a spherical 
geometry which may not accurately represent the Bl-E conditions. (Note that our dust opacity as¬ 
sumptions produce more significant uncertainties than adopting the dust temperature as a proxy of 
the kinetic gas temperature; see Appendix|Al) Nevertheless, B1-E2 has an enriched chemistry with 


highe r abundances of nitrogen-bearing molecules and carbon chain molecules (see iSadavov et al 
2013 ), where the former are generally associated with the freeze-out of and its isotopo- 


logues, and the la tter are generally thought to form on the surfaces of dust grains with carbon-ric h 
icy mantles (e.g., iDi Francesco et al.l 120071 : iBergin Sz Tafallal 120071 : iHerbst Sz van DishoeckIl2009ll . 
though recent studi es have suggested th at gas-phase productions of carbon-chain molecules are 


non-negligible le.g.. iBalucani et al.ll2015lj . 


Based on the combination of lower abundances and more pronounced nitrogen-rich chem¬ 
istry, B1-E2 appears to be uniquely depleted in gas by a factor of two. Depletion processes, 

such as molecular freeze-out, will produce thicker ice mantles on the dust grains toward B1-E2 
comp ared to the other substructu res. Thicker ice mantles should increase the B1-E2 dust opacities 
(e.g., Ossenkopf fc HenninsJ 1994 1. which will consequently lower the measured column densities 
and bring Bl-E2’s X(C^®0) values into agreement with those of the other substructures. Neverthe¬ 
less, B1-E2 has a very moderate decrease in N(C^®0) - only a factor of two - and it isn’t clear how 
much the ice mantles (and thus, the opacity) will be affected from such a small level of depletion. 
Moreover, to have higher opacities from thicker ice mantles, one must still assume that B1-E2 has 
higher levels of freeze-out compared to the other substructures. Therefore, we believe B1-E2 
is depleted in gas by at most a factor of two due to molecular freeze-out. 

The onset of depletion toward a prestellar core precursor has not been previously identi¬ 
fied. This depletion (a factor of two) is much weaker than what is typically dete cted toward dense 


prestellar cores fe.g.. ICaselli et al.lll999l : iBergin et al.l 120021 : iFord fc Shirley 120111 1. and corresponds 


better to the depletion seen in the lo west density objects (ILippok et al 


201 ^ or in very young 


sources (e.g., iTafalla fc Santiago! 120041 ). S ince B1-E2 is slight ly more centrally condensed than the 
other substructures (e.g., see FigurefTland lSadavov et al.ll 2012 l h this distinction in chemistry may be 
tied to its density. Alternatively, B 1 -E 2 is also dynamically distinct from the other substructures, 
suggesting that the evolution of turbulence and chemistry in dense cores appears to he contempora- 
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neous, such that only those objects with low internal turbulence also experience depletions of their 
carbon-bearing molecules and enhancements of their nitrogen-bearing and carbon-chain chemistry. 
We discuss this chemodynamical picture in more detail in Section 15.2.31 


5.2.2. Comparisons with Simulations 


Simulations of core formation have the difficult task of creating small-scale structures within 
larger-scale clouds under a variety of physical conditions (e.g., turbulence, magnetic fields, chem¬ 
ical processes). The high density struct ures that form at early times i n simulations generally 
have supersonic turbu l ent motions (e.g., iBallesteros-Paredes et al.l l2003l: iLi &: Nakamural l2004l : 


Clark &: Bonnelll l2005l : IVazauez-Semadeni et al 


2005 


Ballesteros-Paredes et al 


200711 . Most of 


these early objects are transient and are destroyed within a crossing time unless they are able to 
build up sufficient mass (e.g., through mergers, ambipo lar diffusion) to form long-lasting structures 


( Clark &: Bonnell 


2005 


Vazauez-Semadeni et al.ll2005ll . Although the Bl-E substructures have su¬ 


personic non-thermal line widths as expected in some simulations for high-density objects at early 
times, they also appear pressure-confined (see Section I3.4ji rather than transient. Thus, the Bl-E 
substructures are not at the earliest stages when dense, transient objects hrst condense out of the 
ambient material. Rather, the Bl-E substructures are likely persistent objects, that will eventually 
form into dense (prestellar) cores. 

The high density structures that form out of these simulations can also inter act with each 


other, causing them to evolve dynamically separate from their surrounding cloud. iKlessen et al 


([ 200 ( 11 ) found that mutual gravitational interactions between dense structures could homogenize 
their relative velocities compared to their initial velocity distribution, and that these interactions 
become more effective as the dense structures gain more mass. Thus, it is possible that core 
formation only proceeds when the small-scale, high density structures decouple from the cloud. 
Nevertheless, the expected timescale for these interactions is very long, ~ Assuming a typical 
density of 10^ cm“^ for the Bl-E clump, the expected free fall time is t ff = \/3tt/32Gp ^ 1 Myr, 
and a timescale of ~ is older than all the YSOs in Perseus (e.g., iBallv et al.ll2008l b If core 


formation proceeds only after the early precursors decouple from t he ambient cloud via mutual 
interactions, this process must occur more efficiently. Alternatively, iBailev &: Basul (j2012l j found 
that material with different length scales (e.g., clumps and cores) form separately over different 
timescales. Thus, the Bl-E substructures may have started decoupled from their parent clump, 
such that they formed with more similar relative motions than the random motions within the Bl- 
E. Whichever mechanism decouples high density structures from their parent cloud (e.g., mutual 
interactions or in-situ), it must occur early in the formation of these objects. 

Some simulations have suggested that high density structures inherit th eir kinematic properties 
from the material out of which they fragmented (e.g., iKlessen et al.l 120001 1. In particular for the 
Bl-E substructures, if we assume that their non-thermal NH 3 (1,1) line widths are due to turbulent 
motions (e.g., not multiple subcomponents; see Section [5.2.ip . then their dense gas dynamics, as 
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traced by the NH 3 (1,1) data, have very similar turbulent velocity dispersions as the bulk gas. 
While the Bl-E substructures appear decoupled from their parent clump, their internal kinematics 
still resemble that of Bl-E, suggesting that the timescale associated with a core losing its internal 
turbulence may be much longer than the timescale associated with decoupling dense structures 
from the ambient cloud. 


The substructures also likely inherited chemical properties from the Bl-E clump. Most of the 
Bl-E substructures are chemically unevolved compared to dense (prestellar) core s in that they show 
neithe r or depletion nor enriched nitrogen or carbon-chain chemistry. iTafalla &: Santiago 

(|2004ll attributed low CO depletion to very young objects, such that proc esses like ambipo l ar dif - 
fusion have insufficient time to affect the internal chemistry. Similarly, iNakamura fc: Lil (j2005l l 
attributed low CO depletion levels to formation mechanisms like rapid turbulent compression, 
where an object forms too quickly for its carbon-bearing molecules to sufficiently deplete. Al¬ 
though, we find no in dication of C^^O dep letion in most s ubstructures, we do not see a rich carbon 
chemistry (e.g., CCS dSadavov et al.ll2ni2l ) as predicted bv INakamura fc Lil (j2005|). Thus, the Bl-E 
substructures appear to be very young. 

On the other hand, B1-E2 is the only source to show signs of slight C^®0 depletion and mod¬ 
erate enriched nitrogen-chemistry. Since B1-E2 is also the least turbulent substructure, chemical 
processes may also depend on the local turbulence. Previous studies have connected turbulence 
and cloud chemist ry on scales of the mixing length, though the spatial extent of turbulent mixing is 


still debated fe.g..lCanto fc Raga 


2003; 


Martinell et al 


boofil l. 


1991 


Xie et al. 


1995; 


Bawlings fc Hartquistlll997l : lYate fc Millar 


Most of these studies, however, examined turbulent transport across 
clouds, whereas for Bl-E, we need to determine whether or not turbulent mixing can replenish 
the carbon-chemistry over the scales of a cold, supersonic “core.” In addition to turbulent mix¬ 
ing, s hocks can replenish carbon-bearing molecules from icy grain mantle s back into the gas phase 
(e.g., IPraine et al.lll983l; IWilliams fc Hartauistlll984l ; iGuillet et al.l 1201111 . At the expected shock 
veloc ities of 2 — 3 km s~ ^ ( Pon et al.ll2014l ). dust grain mantles are unlikely to be entirely removed 
(e.g.. iDraine et al.lll983l ). Instead, gas-grain interactions involving h eavier molecular gases like CO 
can sputter dust grains or their icy mantles at these lower velocities (j.Iimenez-Serra et al.ll2008l ). In 
the sputtering case, a smaller fraction of gas-grain interactions will result in CO freeze-out. More¬ 
over, sputtering of dust grains with icy mantles may release additional CO molecules (e.g., those 
molecules that did freeze-out) back into the gas phase, thereby keeping a quasi-equilibrium gas 
phase abundance. Thus, the supersonic Bl-E substructures may have higher C^^O abundances due 
to dust sputtering interactions, whereas the transonic B1-E2 source does not have the kinematics 
for sputtering, and as such, its dust-grain interactions primarily result in freeze-out. Both chemi¬ 
cal mixing and small-scale shock replenishment are possible explanations for the chemistry in the 
Bl-E sources. Detailed chemical models at both clump- and core-scales in a supersonic turbulent 
environment with low-velocity, small-scale shocks are needed to study the adsorption of CO ices 
and the transport of carbon-rich material to core scales. 
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5.2.3. Connecting Dynamics and Chemistry 

Our observations suggest the following chemodynamical picture for the formation of relatively 
isolated prestellar cores; (1) core precursors condense out of their turbulent parent clump with their 
kinematic and chemical properties inherited from that medium. These core precursors either form 
decoupled from their surroundings or are promptly decoupled, e.g., via efficient mutual gravitational 
interactions; (2) the highly turbulent core precursors dissipate their turbulence over some time, e.g., 
via small-scale shocks, while retaining their carbon-rich chemistry; (3) as the turbulence dissipates 
and the density increases, carbon-bearing molecules become depleted from the gas phase enriching 
nitrogen-bearing molecules and large carbon-chain molecules; (4) the core precursors continue to 
evolve on this path of turbulence dissipation, carbon depletion, and chemical enrichment, eventually 
forming a quiescent, chemically layered prestellar core. Whether the turbulence dissipation and 
chemical depletion/enrichment occur simultaneously or after a threshold low turbulence level is 
reached is still unclear. Thus, most of the Bl-E substructures are in the above Stage 2, since they 
have not yet dissipated enough turbulence to allow for significant depletion (i.e., they are 

still carbon-rich). The B1-E2 substructure, however, is slightly more evolved in Stage 3. The C^®0 
depletion and chemical enrichment detected toward this object are both weak compared to most 
prestellar cores, indicating that B1-E2 is still very young. Thus, the Bl-E substructures represent 
excellent candidates for core precursors and a very good opportunity to catch the initial stages of 
core formation. 


6. Conclusions 


We have presented new and observations in the J = 1 — 0, J = 2 — 1, and J = 3 — 2 
transitions of the Perseus Bl-E clump. To compare all data equally with the Herschel-deiwed dust 
temperatures and column densities, we convolved the molecular line observations to a common 
resolution of 36.3'h We compared and contrasted the ki nematics of th e entir e Bl-E clump with 
those of the individual nine substructures identihed by ISadavov et al.l ([2013) • Additionally, we 
characterized the and chemical signatures. Our main conclusions are: 


1. The Bl-E clump has a velocity gradient of 1 km s ^ pc ^ detected in both the (1-0) 

and C^®0 (1-0) spectra. This velocity gradient indicates that the gas velocity is increasing 
from the NW to the SE at an angle of ~ 20° — 30°, north to west. In contrast, material in 
the larger Perseus cloud has velocities generally increasing from the SW to the NE. Since 
differences in gas ffows can lead to cloud interactions, the formation of substructures may 
be triggered in Bl-E by the shear of these different ffows or by compression if the ffows are 
colliding. At this time, the origin of this different velocity gradient and its impact on Bl-E 
is unclear. 
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2. We find that the Bl-E substructures have more homogenized motions than typical turbulent 
velocities detected in the larger-scale ambient c lump (as traced by C^^O), similar to what has 
been found in other clumps in Perseus (e.g., see lKirk et alJl2009l : iKirk et alJl201Cll ). Moreover, 
the NH3 emission does not show a velocity gradient. Thus, the Bl-E substructures appear 
decoupled from their environment, indicating that dense cores either form decoupled from 
their parent clump or obtain such homogenized motions through external factors early in 
their evolution. 


3. We find that the turbulence in Bl-E is difficult to quantify from the and spectra 

as both may contain multiple velocity components. Using single Gaussian fits, the median 
non-thermal velocity dispersion of the ( 1-0) spectra is ^ 0.4 km s“^, which agrees 

well with the broadest NH3 (1,1) emission from ISadavov et al.l (j2012l ). We propose that core 
precursors form with the same turbulent properties as their parent clump. As thes e sources 
evolve, their turbulence dissipates (e.g., via small-scale shocks; see lPon et al.ll2014l ) and the 
core precursors become quiescent. Nevertheless, high resolution observations are necessary 
to determine whether the broad line emission toward the Bl-E substructures are due to 
supersonic turbulence or multiple subcomponents along the line of sight. 


4. Using the radiation transfer code, RADEX, we find that Bl-E is best fit with X(C^®0) 


2 X 10 For N(H2) < 8 x 10^^ cm we see primarily lower values of X(C^^O), suggesting 

TCI cn rvn T "fi ttP (tI- -f-V»/-»ci/-» i n t-v-i /-IjtitI- i/-vo A "P 1 


that selective dissociation of is significant at these column densities. At N(H2) > 10 

cm“^, we find that the column density toward the B1-E2 substructure is uniquely lower 
than the other sources. We consider B1-E2 to be depleted in by a facto r of two, in 

agreement with its relatively enriched chemistry, as seen bv ISadavov et al.l (|2012l ). 


5. Most Bl-E substructures are undepleted in and without enriched chemistry. The excep¬ 
tion, B1-E2, is also the only substructure with quiescent (transonic) turbulence. Therefore, 
we suggest that the depletion of carbon-bearing molecules and the enrichment of nitrogen¬ 
bearing molecules are concurrent with turbulence dissipation, where cores themselves evolve 
from highly turbulent, undepleted objects to less turbulent objects with CO-depleted chem¬ 
istry. 


The Bl-E clump is an excellent laboratory to explore core formation. Since this region has 
not formed a previous generation of cores or protostars, we have the unique opportunity in Bl-E to 
observe the initial conditions associated with dense cores condensing out of their parent clump. In 
particular, our observations have hinted at a number of interesting possibilities, including possible 
triggering (from shearing or colliding flows) and an evolutionary sequence connecting turbulent 
dissipation to chemical enrichment in the core precursors. Nevertheless, further study of this region 
is necessary to understand better the processes influencing this clump. In particular, simulations 
that explore the core formation process and chemical depletion in a highly turbulent environment 
are key to connecting the initial characteristics of these objects to the characteristics of their 
surrounding clump environment. Additionally, studies that explore the impact of low-velocity cloud 
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interactions are necessary to determine the significance of a clump with a velocity flow different 
from that of its parent cloud. With such studies, we may have, for the hrst time, the necessary 
tools to constrain core formation. 
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A. Justification of Dust Temperature 


In Section [H we used the radiative transfer code, RADEX, to model the multi-J ^^GO and 
G^^O observations using the Herschel-deiived dust temperature as a proxy for the kinetic gas 
temperature. Since the and line emission mainly trace material at n < 10^ cm“^, the gas 
and dust may be decoupled such that the kinetic gas te mperatures (at the densities traced by ^^GO 


and C^®0) are warmer than the dust temperatures le.g.. iGoldreich fc: Kwanlll974l : lGoldsmithll2001 


This preprint was prepared with the AAS lAERK macros v5.2. 
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Di Francesco et alJ 120071 : ICeccarelli et alJ 120071 ). Morevoer, the best-fit radiative transfer models 
had different excitation temperatures for each J-line, suggesting the molecules whose emission we 
observe are not in thermal equilibrium with the unknown kinetic gas temperature. 


We estimated the kinetic gas temperature using the COMPLETE survey (jBidge et akl 120061 ). 
which mapped th e entire Perseus mole cular cloud in (1-0) and (1-0) at ~ 45'' resolution. 

Using these data. iPineda et al.l (2003) determined the excitation temperature, T^x, from the 
(1-0) data assuming optically thick emission, following; 


ATmb = (1 - e-")^ [f{%x) - f{TcMB )], 
k 


(Al) 


where ATmb is the main beam brightness of the line, r is the optical depth, h is the Planck 
constant, k, is the Boltzmann constant, and f{T) corresponds to. 


fiT) = 


exp 


hiy 

kT 


- 1 


-1 


(A2) 


and f{T) is measured for the excitation temperature {T^x) and the background temperature, given 
by the cosmic microwave background (CMB). Assuming a CMB temperature of 2.73 K, Tg,x is 
measured as, 

(A3) 


T — 

-‘-ex — 


In 1-1- 


5.53 K 

ATmb + 0.84 K 


Pineda et al.l (j2008li reported typical Tf,x values of T^x ~ 11 K in different regions of Perseus at 
5' resolution (to match the resolution of their extinction data). We followed the same procedure to 
characterize T^x on a pixel-by-pixel basis at the native resolution of the COMPLETE data (~ 45"). 
We used single Gaussian fits to measure the peak line brightness (Tmb) at each pixel, assuming a 
beam correction factor of ~ 0.45 to convert from TJ to Tmb- We also assumed the emission filling 
fraction is unity. 


Figure [A. II shows the excitation temperature across Bl-E determined from Equation IA3I with 
the integrated intensity map of the ^^CO (1-0) observations from COMPLETE. For comparison, we 
also show the Hersc hel-deiived dust temperature {T^ust) map from SED fitting (see Section [2] and 
Sadavov et al.ll2014l for details) and the ratio of T^ust and T^x- Both temperature maps show very 
different structures. The ^^CO-derived T^x measurements range over 10 — 22 K, with the highest 
values towards the center of the field, and no correlation is found between Tga, and column density 
(contours). Since the ^^CO (1-0) line emission is optically thick, it will not trace the dust emission 
(e.g., see the ^^CO integrated intensity map in Figure IAT]) . In contrast, the iLersc/ief-derived dust 
temperatures range over 12 — 16 K, with a strong correlation between T^ust and column density 
such that the lowest temperature values are towards the densest regions of Bl-E. Similarly, the 
^^CO and C^^O integrated intensity maps also trace the dust emission well (see Figures [3] and |1|). 
Therefore, T^ust is likely a better representation of the kinetic gas temperature associated with the 


i3r 


CO and C^®0 gases than the 


12 


CO-derived T^x values. 
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Fig. A.l.— Top: Excitation temperature (Tex) toward Bl-E (left) and integrated intensity of 
(1-0) li ne emission (right) . Both maps are at 46 ” resolution and used data from the COMPLETE 
survey (jRidge et al.l 1200611 . Bottom: Herschel-deTwed dust temperature {T^ust) at 46" resolution 
(right) and the ratio of T^ust and T^x with column density (left). We used modihed blackbody 
SED fitting to determine T^ust from 160 — 500 fim Herschel data assuming a dust emissivity index 
of /3 = 2. The T^ust map has also been convolved to 46" resolution to match the resolution of 
the (1-0) data. Contours show Herschel-derived column density levels from HGBS data of 

5.0 X 10^^ cm“^,7.0 x 10^^ cm“^,8.5 x 10^^ cm“^,and 10.5 x 10^^ cm“^. 


Although the gas and dust are likely not coupled toward most of Bl-E, the kinetic gas temper¬ 
ature and t he dust temperatur e are expected to vary by only a few Kelvin at densities of n 103-4 
(e.g., lYoung et al.ll2004l : ICeccarelli et al.l 120071 ). which represent the bulk densities of Bl-E 
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fsee ISadavov et alJl2012l ). Therefore, we re-ran the RADEX models for a randomly selected sample 
of test pixels assuming the CO gas kinetic temperature is warmer than the dust temperature by 
2 K. We found that the column densities were generally lower by < 20%, whereas the 

column densities were lower by < 10%. These uncertainties are relatively minor, considering the 
complex spectra of the and gases and the simple geometric assumptions in the RADEX 
routine. As such, in the absence of proper gas kinetic temperatures, the Rersc/ieZ-derived dust 
temperatures give a good first-order estimate for abundances from RADEX. 





